Realistic shell model calculation of 2v/3/3 

nuclear matrix elements 
and role of shell structure in intermediate 

states 

H. Nakada 

Department of Physics, Faculty of Science, Chiba University, 
Yayoi-cho 1-33, Inage-ku, Chiba 263, Japan 

T. Sebe 

Department of Applied Physics, College of Engineering, Hosei University, 
Kajino-cho 3-7-2, Koganei, Tokyo 184, Japan 

and 
K. Muto 

Department of Physics, Faculty of Science, Tokyo Institute of Technology, 
Oh-okayama 2-12-1, Meguro-ku, Tokyo 152, Japan 



PACS numbers: 23.40.-s, 21.60.Cs, 23.40.Hc 



Keywords: 2v(3(3 decay, nuclear matrix element, shell model calculation, shell 
structure 



Abstract 

We discuss two conditions needed for correct computation of 2v(3(3 
nuclear matrix-elements within the realistic shell-model framework. 
An algorithm in which intermediate states are treated based on White- 
head's moment method is inspected, by taking examples of the dou- 
ble GT+ transitions 36 Ar -► 36 S, 54 Fe -» 54 Cr and 58 Ni -> 58 Fe. 
This algorithm yields rapid convergence on the 2v(3(3 matrix-elements, 
even when neither relevant GT + nor GT~ strength distribution is 
convergent. A significant role of the shell structure is pointed out, 
which makes the 2v[5[5 matrix-elements highly dominated by the low- 
lying intermediate states. Experimental information of the low- lying 
GT ± strengths is strongly desired. Half-lives of T^EC/EC; 36 Ar -> 
36 S) = 1.7 x 10 29 yr, Tfj 2 (EC /EC; 54 Fe -► 54 Cr) = 1.5 x 10 27 yr, 
Ty 2 (EC /EC; 58 Ni -► 58 Fe) = 6.1 x 10 24 yr and T?J 2 (P+ /EC; 58 Ni -» 
58 Fe) = 8.6 x 10 25 yr are obtained from the present realistic shell-model 
calculation of the nuclear matrix-elements. 



1 Introduction 



The nuclear double-/? ((3(3) decay is a good probe to investigate neutrino 
masses and some basic properties of the weak interaction |l|, |2|, |3|. The Ou 
mode, if it exists, immediately indicates a defect of the standard model. 
On the other hand, the 2v mode occurs within the standard model. This 
mode has a practical importance, because it governs the lifetime of the nuclei 
which concern (3(3 decay. Moreover, the observation of the 2v(3(3 decays|| has 
revealed serious discrepancies between earlier predictions and the observed 
half-lives. It probably originated in a problem of the nuclear structure theo- 
ries. This problem may be crucial also to the 0v(3(3 decays, because some of 
the nuclear matrix-elements associated with the Ou mode have similarity to 
the 2v mode. 

The 2v(3(3 decay takes place by two sequential Gamow- Teller (GT) tran- 
sitions, via virtual intermediate states. It is not always possible for the cur- 
rent nuclear structure theories to give reliable predictions even on single-GT 
transitions, since the GT transitions are sensitive to some details of nuclear 
many-body wavef unctions. The most promising approach to the GT tran- 
sitions seems to be realistic shell-model calculations. It is yet difficult to 
apply realistic shell-model approaches to the relatively heavy nuclei where 
the 2v(3(3 decays have been observed. It may become possible in the near 
future, however, to carry out a reliable shell-model calculation in a few of 
these nuclei, with assistance of the growing computer power. 

Several realistic shell- model calculations have been reported for the 48 Ca — > 



Is 



Ti (3 (3 decayP, 10, 11], although it has not been observed. In Ref. [ 
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it has been pointed out that, among a large number of possible intermedi- 
ate states, the lowest 1 + state gives a dominant contribution to the 2vf3f3 



process. In Ref.[10|, a method to handle the intermediate states has been 



developed based on Whitehead's moment method|12[. Rapid convergence by 
this method has been shown for the 48 Ca decay. The 48 Ca case is, however, 
somewhat special, as will be discussed later. It is worth examining generality 
of these suggestions, by taking other sample nuclei. 

While most experimental efforts have been cast on the double-/? - ((3~ (3~) 
decays, the (3 + side is also important. Since the electron capture (EC) and 
positron emission (f3 + ) are possible in a single-GT + transition, we have three 
cases in the double-GT+ transitions; EC/EC, (3 + /EC and §. It 

has been pointed out that, in the Oz/ mode, they seem to yield quite a 
different constraint on the neutrino masses and right-handed weak-current 
parameters ||, although the double-GT + transitions are more difficult to be 
observed than the (3~ (3~ cases. From the viewpoint of the nuclear structure 
theories, there is no essential distinction in the algorithm to compute the f3(3 
matrix-elements, between the f3~ and the f3 + cases. In the A <^ 60 region, 
where realistic shell-model wavefunctions are available, the (3 + side delivers 
a few candidates for the [3(3 decay in addition to 48 Ca. These candidates 
give us a good opportunity to test the theoretical methods. One important 
issue is practical treatment of the intermediate states. Primarily focusing on 
this point, we shall discuss in this article how to compute the 2u(3(3 nuclear 
matrix-elements within the realistic shell-model framework, by taking a few 
sample cases from the f3 + side; the 36 Ar — > 36 S, 54 Fe — > 54 Cr and 58 Ni — > 58 Fe 
decays, where realistic shell-model wavefunctions are available for the initial 



and final states]?], |J. We will also discuss the key features to describing the 
2vf3(3 nuclear matrix-elements. 



2 Shell model calculation of 21/(3(3 matrix ele- 
ments 

In this section, we present a realistic shell-model calculation of the 2vj3j3 
nuclear matrix-elements. We discuss the conditions required for treating the 



intermediate states. The algorithm proposed in Ref . [|10j is argued in the light 
of these conditions, and is re-examined for the examples of the 36 Ar — > 36 S, 
54 Fe _^ 54 Cr and 58 M _^ 58 Fe decays. 

With respect to the 36 Ar decay, a realistic shell-model hamiltonian has 
been established in the sd-shell region; the so-called USD interaction |7J. The 
54 Fe and 58 Ni nuclei belong to the middle pf-shell region. Though diago- 
nalization of a hamiltonian in the full pf-shell is not yet possible for the 
middle pf-shell nuclei, a large-scale calculation in a loosely truncated space 
is successful to describe spectroscopic properties||. The k < 2 configu- 
ration space, where k represents number of nucleons excited from O/7/2 to 
(O/5/2 IP3/2IP1/2), is adopted in this calculation, together with the Kuo-Brown 
hamiltonian]!^]. The 2vj3f3 nuclear matrix-elements contain a product of a 
GT + matrix-element from the initial state and a GT~ matrix-element from 
the final state. In the 54 Fe and 58 Ni decays, the k > 2 configurations are 
forbidden in the intermediate GT + statefll^l, as far as the initial state is 



restricted to the k < 2 configurations. Notice that the lowest (i.e., k = 0) 
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configuration shifts from the parent nucleus to the intermediate one, because 
the maximum number of nucleons which can occupy the O/7/2 orbit changes 
from nucleus to nucleus in this region. Unless a configuration has both the 
GT strengths from the initial and final states, this configuration does not 
contribute to the 2u/3j3 matrix-elements. Accordingly, only the k < 2 con- 
figurations are relevant to the 2v(3(3 process, also in the GT~ strength from 
the final state. This situation is favorable because the k < 2 wavefunctions 
have been inspected well by spectroscopic studies, including some single-GT 
transitions [Q. It should be stated here that these shell-model calculations 
in the sd- and pf-shell do not provide us with correct total binding ener- 
gies. The Coulomb energy is neglected in the sd-shell hamiltonian, while the 
ground-state energy has been out of scope in the pf-shell calculation. We 
shift the energies of the intermediate states so as for the ground-state energy 
of the intermediate nuclei to fit to the experimental value, relative to the 
initial or final states. In the 36 Ar — > 36 S and 54 Fe — > 54 Cr decays the modes 
other than the EC/EC mode cannot occur because of the Q- values, while 
the /3 + /EC mode has positive Q-value for 58 Ni — > 58 Fe. 

Though the practical application will be devoted to the f3 + side, we here 
discuss how to calculate the 2vf3f3 nuclear matrix-elements more generally. 
The nuclear matrix-element of the 2vf3f3 process is written as 

M„(2^)^ </ " r(GT ; ) ""' ) <1 r(GT±) " i) . (U 

where Ei and E m denotes energies of the initial and intermediate nuclear 
states, respectively. The symbol to expresses energy of leptons emitted (or 
captured) during the step producing the intermediate state from the initial 
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state. Its value ranges from zero to the Q-value of the (3(3 process. The 
initial state (|z)) is the ground state of the parent nucleus. In order for 
the (3(3 decays to be detectable, a single-GT decay from the parent nucleus 
must be forbidden. This happens for some decays from even-even nuclei, 
owing to the pairing correlation. Hence the initial states of the (3(3 decays 
have J p = + . The intermediate nuclear states {\m)), whose eigenenergy is 
denoted by E m , should have J p = 1 + . Although the final state (|/)) may 
have J p = 2 + in some (3(3 decays, the ground + state is expected to give the 
dominant branch IJ. In the practical applications to be discussed, the final 
states must be the ground state, because of the Q-values. The contribution of 
the Fermi-decay mode due to the isospin-mixing is known to be negligible^], 
and is omitted in Eq. ([]]). The bare GT operator is given by 

I *»(GT ± )=X;^ < t tl4> (2) 

i 

where the summation runs over constituent nucleons and g A = Ga/Gv = 
1.26. It has been known that, however, effects of the core-polarization (CP) 
and meson-exchange currents (MEC) should be taken into account in shell- 
model calculations on the GT processes [|l^, For this reason we use effec- 
tive GT operators, 

T(GT ± ) = T free (GT ± ) + 5T(GT ± ), (3) 

where 

5T(GT±) = £ {5g A (nl)<Ti + Sg lA (nl)k + 5g pA {nl ) n'l')[Y i - 2 \h0a l }^}t ± ,. 

i 

(4) 



In the shell-model approach, the summation in Eqs. (0) and (^) is restricted to 
valence nucleons. The parameters Sg^, SgiA and 8g p A depend on n and I; the 
quantum numbers of single-particle orbit which the i-th nucleon occupies. In 
the sd-shell region, an effective operator for the USD wavefunctions has been 
obtained] 17| , by making a global fitting of the single-particle parameters of 



Eq.(|4]) to the measured -B(GT) values. In the pf-shell region, we adopt the 
parameter-set evaluated by Towner from microscopic standpoints |T3|, which 



has been shown to yield a sound agreement with the data|T4|, combined with 
the k < 2 Kuo-Brown wavefunctions. 

Application of Eq. QT]) associated with the above GT operators to the 2v/3/3 
matrix-elements implicates the following two assumptions: (i) The contribu- 
tion of high-lying intermediate states beyond the Ohoo space is negligible, 
although a certain amount of the GT strength is carried out of the Ohu 
space due to the CP mechanism. This is expected because, in addition to 
large energy-denominator, the high-lying GT strengths from the initial and 
final states are unlikely to be coherent, and has been supported by a statisti- 
cal estimate [18|]. (ii) The two sequential GT transitions are separated so well 



that the intermediate nuclear states should be the nuclear energy-eigenstates. 
The time interval of the sequential GT transitions is ruled by the energy de- 
nominator via the uncertainty principle, leading to (1 — lOMeV) -1 . It is much 
longer than the time-scale of the weak current, which is the inverse of the 
gauge-boson mass, ~(10 2 GeV) _1 . The MEC process takes (10 2 MeV) _1 at 
most, owing to the pion mass. Because the MEC concern the <5T(GT) term, 
the correction due to the time of the meson propagation will not exceed a 
few percent of the total 2h>/3/3 matrix-elements. Contribution of other higher- 
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order diagrams can be estimated to be negligibly small, in a similar manner. 
The assumption (ii) is thus plausible within a few percent accuracy. 

The explicit construction of the intermediate energy-eigenstates (\m)) in 
Eq. ([I]) is difficult in most cases, even if the wavefunctions of the initial and 
final states are obtained. Although the method of solving a coupled linear 



equation JTTj] is less elaborate than full diagonalization of a hamiltonian, still 
it is not easy to obtain all the possible intermediate states. This is because 
there may be quite a large number of intermediate states. In the closure 
approximation adopted in earlier studies||, E m of Eq. ([!]) is replaced by a 
constant average value. We then have, assuming a + final state, 

' -db\ TVOrpi 



= ~B — {f\\T{GT ) ■ T(GT ± )\\i), (5) 

and do not have to handle the intermediate states explicitly. The average 
energy has been estimated mainly from the GT^ 1 resonance energies. It 
has been recognized, both theoretically and experimentally, that this closure 
approach hardly gives reliable matrix elements. This indicates a problem in 
the estimate of E av and/or in the procedure taking energy average itself. 

However, averaging of energy is justified if we carry it out in a sufficiently 
small energy range. Let us consider a set of 1 + states of the intermediate 
nucleus {\n}} which exhausts the relevant GT^ strength from the initial 
state. The states |n)'s do not have to be energy eigenstates. It is possible to 
insert Y^n |™)(™| into the numerator of Eq. ([[]), 

m.n E m — Ei + U) 



It should be commented that the set {|n)} does not have to run out all the 
GT strength; it has only to exhaust the relevant strength. For instance, 
even if the GT strength has a certain isospin distribution, only the lowest 
isospin component needs to be considered for because higher isospin 

components do not have GT-transition strength to |/). A state \h) has large 
overlaps with some eigenstates and negligibly small ones with others. These 
overlaps depend on E m . It is postulated that the energy distribution of |n) is 
peaked so sharply that |m)'s with large (m\n) could have eigenenergies close 
to E n = (n\H\h). We then obtain, via the approximation E m ~ E n under 
the presence of (m\n), 

Mu&ffl * £ p \r— E(/l l T ( GT± ) H m ) Hn) (*\\T (GT ± ) 1 1<>. (7) 

The closure with respect to \m) is available at this stage, giving 

n E n — Ei + U) 

The state \n) plays a role of a doorway state for a certain number of energy 
eigenstates \m). In comparison with Eq. ([!]), Eq.(|) appears simply an re- 
placement of \m) by \h) and E m by E n . Namely, a good approximation will 
be obtained by using an appropriate set of doorway states instead of the in- 
termediate energy-eigenstates. The approximation of Eq.(|]) can be regarded 
as an extension of the closure approximation. It indeed reduces to a closure 
approximation when a single intermediate state (the GT state) is used for 
|n). The remaining problem is how to generate {\n)} efficiently 

In Ref.fU, a cancellation mechanism has been shown by using a weakly- 
coupled two-configuration model. It is noted that this mechanism is compre- 
hended in the above equations; the essential part of the argument in Ref.[H 



reappears if we regard \n) as a state having a considerable GT 1 * 1 strength 
from \i) but no GT T strength from |/). The vanishing GT T strength pro- 
hibits \n) from contributing to M tJ (2uPP) (see Eq.(|8|)). There may be one 
or more states with opposite character, having a GT T strength from |/) but 
no GT 1 * 1 strength from Suppose that some of these states have energies 
close to Ek, and couple weakly to \n). This coupling gives rise to a certain 
fragmentation of |n) over several eigenstates, which is here denoted by \mi), 
\17i2), ■ ■ •, \m r ). Their eigenenergies are close to E n . Each of \rrii), \rri2), ■ ■ ■, 
\m r ) may have a sizable contribution to M UJ {2vf3j3) 1 because they have GT 
strengths both from \i) and |/). These contributions, however, cancel one 
another, and the result obtained solely by \n), which gives the vanishing con- 
tribution to M LJ (2i'P/3), is recovered. The whole contribution of \mi), \rri2), 
■ ■ -, \m r ) to M w (2vfiP) is represented by \n), as is shown from Eq.(||) to (||) 
in more general cases. 

The set {\n)} has been assumed to exhaust the transition strength from 
the state \i), in the above discussion leading to Eq.fl8|). It should be noted 
that an analogous equation is derived by adopting a set of doorway states so 
as to exhaust the relevant GT strengths from |/). These two ways give the 
same value if they are convergent, although the speed of convergence may be 
different. 

The present derivation of Eq. (|8D discloses what conditions are required for 
{|fi)}. There are two conditions: (a) The set {|n)} must exhaust the relevant 
GT strength either from the initial or final state, and (b) each doorway state 
\n) has an energy distribution over \m) with sharp peak around E%, so as 
for E m ~ Ea in the energy denominator to be justified. A promising way to 
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obtain such approximate intermediate states has been proposed in Ref . [TO] , 



by applying Whitehead's moment method|L2|]. The algorithm is as follows: 



(I) Produce the state exhausting the relevant GT 1 * 1 strength from the initial 
(final) state, |GT^) oc P rel • T(GT ± )|i(/)), where P re i stands for a 
projection operator which picks up the relevant configurations. 



(II) Generate a set of bases by operating the shell-model hamiltonian H on 



(III) Diagonalize H within the space = {|GTtJ /} ), H\GT± {f) ), • • • , H N \GT± {f 
Each energy eigenvalue and eigenstate in this subspace are taken as 

and \h). 

(IV) Increase N, until 71^,(2^/3/9) converges. 

Although the steps (|I|) and (PTl] ) are common with the Lanczos diagonal- 
ization procedure, the convergence of M UJ {2vj3j3) does not necessarily imply 
the convergence about the intermediate states, as will be shown later. The 
condition (a) is obviously fulfilled, since contains jGT^), into which 
all the relevant GT strength is concentrated. The simplest case of N = 1 is 
equivalent to the closure approximation of Ref. Speed of the convergence 
for N depends on how well (b) is satisfied. This point will be examined below 
for the practical cases. 

Since all the (3(3 decays to be investigated belong to the j3 + side, the 
present practical calculation has to do only with T(GT + ), not T(GT~), in 
Eqs.(|TJ — |S|). For the GT state generated from the initial or final state at the 
step (|), |GT+) or |GTj) is considered. When |GT^) is chosen at (|), we 
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do not use the P re i operator (i.e., P re \ = 1). For |GT7), P re \ represents the 
isospin projection for the 36 Ar decay, while for the 54 Fe and 58 Ni decays the 
projection onto the k < 2 space is included as well. 

Because of the phase space of the emitted leptons, uj ~ Q/2 yields a 
main contribution to the (3(3 probabilities, where Q denotes the Q- value of 
the (3(3 process. We examine convergence of the above method, by taking 
uj = Q/2 as a typical value of uj. Fig.|I| shows the convergence for N (the 
number of the iteration). The result on the EC/EC mode is presented for the 
58 Ni decay. The (3 + /EC mode has no essential difference. The number of the 
J p = 1 + , T = 1 (i.e., relevant) states with the s<i-shell configuration is 54 in 
36 CI, the intermediate product of the 36 Ar decay. Hence a full diagonalization 
of the shell-model hamiltonian is possible. On the other hand, for 54 Mn 
and 58 Co, there are more than 10 4 states having J p = 1 + and the lowest 
isospin, even in the present truncated configuration space. Nevertheless, the 
convergence is so rapid that N = 10 could be almost sufficient for all these 
cases. This situation is quite similar to the 48 Ca decay reported in Ref.[[l(] . 
The cases starting from the GT + state and the GT~ state are compared 
for the 36 Ar decay. We view that the GT + case yields somewhat better 
convergence, though the difference is not quite big. The same holds for the 
54 Fe and 58 Ni decays. The reason for this will be discussed later. 

In Fig|| the GT + strengths from 54 Fe are depicted for N = 10, 30 and 
150. Obviously the strength function does not converge with iV = 10 and 30. 
The convergence of the GT~ strength function is somewhat worse than the 
GT + one. Nevertheless, M u) (2u(3(3) is convergent, as shown in Fig.|I|. This is 
because the condition (b) for the set {\n)} is satisfied well, even with such a 
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Figure 1: 



Figure 2: 



small number for N. We can track how each strength for N = 10 diffuses as 
N increases. For instance, about 90% of the GT + strength at E x = 2.7MeV 
appearing in the N = 10 case remains in E x = 2.0 — 3.1MeV in the N = 150 
result. For the strengths with N = 30, about 90% of the E x = 2.9MeV 
strength is shared only by the states at E x = 2.7 and 3.1MeV of the N = 150 
result. In this manner, the N = 10 or 30 strengths, particularly the low-lying 
ones, are found to stay in a small energy range around the original position. 
It is confirmed that the E x < 5MeV strengths are convergent with iV = 150, 
whereas higher-lying ones are not. 
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3 Dominance of low-lying states in interme- 
diate nucleus 



We next investigate, taking notice of the excitation energies, which interme- 
diate states give significant contributions to the 2vf3f3 matrix-elements. 

Fig|| shows contribution of the intermediate 1 + states to M^ = q/ 2 (2^/5/3), 
as a function of the excitation energy. Summed value of the 2uj3/3 matrix- 
elements up to the excitation energy E x , 

(/||T(GT ± )||^)(fi||T(GT ± )|K) 
^ E~ - E- + uj 

n for E A -E <E X ^ n Lj^u 

is displayed, in order to make the convergence for the energy transparent. 
The ground-state energy of the intermediate nucleus is denoted by Eq here. 
Thus the contribution of each intermediate state is presented by the stepwise 
increase or decrease. The result of the full diagonalization is shown for the 
36 Ar decay, while the N = 150 results for the 54 Fe and 58 Ni decays. The 
computation of M u (2uf30) begins with |GT+) at the step (|) (see Section 0). 
It is found that low-lying intermediate states dominate M w (2^/3/3). The 
E x <> 5MeV intermediate states govern M UJ (2i , f3f3) for the 54 Fe and 58 Ni 
cases, while a few of the 5 < E x < lOMeV states give a sizable contribution 
for 36 Ar. This difference in energy might be caused by the difference in 
mass region. Furthermore, the number of the intermediate states giving 
considerable contribution to M w (2v(3(3) is not very large. The dominance of 
a relatively small number of low-lying states in M u (2v/3/3) seems consistent 
with the suggestion in Ref.|| for the 48 Ca decay, although the lowest 1 + state 
is not enough to approximate M u (2u/30) in the present cases. 
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Figure 3: 



We shall look into the 36 Ar — > 36 S case in more detail. The shell- model 
hamiltonian is fully diagonalized, without any further truncation in the sd- 
shell. The energy distribution of the contribution to the EC/EC nuclear 
matrix-element has been shown in Fig|| being repeated in the middle sector 
of Fig.[|. The upper (lower) sector of Fig.f| displays the energy distribution 
of the GT + (GT~ with T = 1) strength from the initial (final) state. The 
GT + strength distribution correlates rather well to the steps of increase or 
decrease in the graph of M [i ,(2i/ / 9/3), though the size of the steps relative to 
the GT + strengths comes smaller as the energy increases. The distribution of 
the T = 1 GT~ strength is quite different from that of GT + , having greater 
weights in high-energy region. Indeed, the overlap between |GT+) and |GTj) 
is very small (~ 2%). Higher-lying GT~ strengths hardly contribute to 
M UJ {2vf3f3) 1 on account of the small overlap to the GT + strength as well as of 
the large energy-denominator. Only a small portion (low-lying part) of the 
GT~ strength contributes to the 2vf3f3 matrix-element. This fact accounts 
for the slower convergence from |GT7) than from |GT+). This mechanism 
applies also to the 54 Fe and 58 Ni decays. 
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Figure 4: 



Contribution of the low-lying intermediate states is enhanced to a certain 
extent, due to the energy denominator. The effect of the energy denomina- 
tor is, however, insufficient to account for the remarkable dominance of the 
low-lying states clarified presently. The particular significance of the low- 
lying GT strengths is qualitatively understood from a general argument on 
the shell structure of nuclei. To the first approximation, the initial and fi- 
nal states have configurations depicted by Fig.|-a) and b), because they are 
lowest-lying. The boxes in the graphs indicate the orbits occupied commonly 
between the initial and final states. The GT + transition from the configu- 
ration of Fig.|5|-a) generates the configurations of Fig.|5|-c), while Fig.|5|-b) 
leads to Fig.||-d). The 2u(3j3 decay proceeds through the common configu- 
ration between Fig||-c) and d). It is immediately found out that only the 
first graphs, which actually imply the low-energy part of the GT + and GT~ 
strengths, are common between c) and d). Therefore, contribution of higher- 
lying GT strengths is hindered by the shell structure, as well as by the energy 
denominator. This shell effect is ignored in the closure approximation, mak- 
ing the energy denominator too large and giving rise to an underestimate of 
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Figure 5: 



M w (2v/3/3), as viewed in the N — 1 limits of Fig.|l| It is also noticed that, 
compared with the GT + transition, the GT~ one has more patterns of high- 
lying configurations, originating in the neutron excess. The central energy of 
the GT~ strength is therefore higher than that of GT + , as has been seen in 
Fig.|], even if we restrict ourselves to the lowest isospin component. 

Let us consider the shell-structure effect on M^{2v(5(5) more concretely, 
for the example of the 36 Ar decay. The lowest configurations of 36 Ar and 36 S 
are 

(p : Od~f 2 , n : 0d~f 2 ) for 36 Ar, (9) 
(p : Orf-/ 4 2 ) for 36 S, (10) 

on top of the 40 Ca core. Thereby the GT + excitation from 36 Ar and the GT~ 
excitation from 36 S primarily produce the following configurations of 36 C1, 

(p : 0d~f 2 , n : 0d~f 2 ), (p : 0d^ 2 0d^ 2 , n : 0d^ 2 ) for GT+ from 36 Ar, 

(11) 

(p : Od"/ 3 , n : 0d~f 2 ), (p : ®d~j 2 , n : 0d^ 2 ) for GT~ from 36 S. 
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(12) 

The common configuration {p : 0d%? 2 , n : 0d^ 2 ), which is the lowest configu- 
ration of 36 C1, gives the main contribution to the 2i>(3f3 matrix-element. The 
(p : 0g^ 3 2 , n : 0d^ 2 ) constitutes only a single 1 + state with T = 1. This com- 
ponent is, however, fragmented in real situations, because other Ohu configu- 
rations appreciably mixes due to the residual interaction. Together with the 
leakage out of the configurations of Eqs.(]9|JT0|) for the initial and final states, 
this fragmentation of the 1 + component distributes the 2uj3f3 matrix-elements 
over a certain range of energy. It is noted that the (p : 0d^ 2 , n : 0d^ 2 ) config- 
uration yields higher energy than the {p : 0d^ 2 0d^ 2 , n : 0d^ 2 ) one, because 
the neutrons occupying Od^/2 are harder to be excited than the protons in 
the same orbit. In terms of the BCS theory, the quasiparticle energy of 0d§/2 
is higher for neutrons than for protons. Hence the GT~ energy should be 
substantially higher than the GT + one. 

In the 54 Fe and 58 Ni cases, the 2uf3(5 processes would be forbidden, if 
merely lj>3/2 were active for neutrons. The principal contribution comes 
from the configurations with the excitation from 1^3/2 to O/5/2. It still holds 
that 2v(3f3 proceeds mainly through lower configurations of the intermediate 
nucleus. Thus the low-lying intermediate states dominate the 2vf3f3 matrix- 
elements, as shown in Fig.|3]. 

We here mention a speciality of the 48 Ca — > 48 Ti decay. To the first 
approximation, 48 Ca can be regarded as a doubly magic core. The lowest 
configuration of 48 Ti is (p : Of 7 j 2 , n : 0fy 2 )- The 2uf3f3 decay is mediated by 
the (p : O/yL, n : 0f 7 } 2 ) configuration, which is the lowest configuration of 
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48 Sc. This configuration allows a single 1 + state, like the 36 Ar case. On the 
other hand, the O/7/2 orbit is isolated and the residual interaction can hardly 
mix this configuration with other ones. Therefore it remains almost pure, 
constituting the lowest 1 + state. This consideration accounts for the reason 
that the lowest 1 + state already gives a good approximation to the total 2vf3f3 
matrix-element in the 48 Ca decay, as pointed out in Ref. ||, unlike the other 
cases investigated in this article. The dominance of the lowest 1 + state has 
been suggested from experiments for the 100 Mo decay [|19|, and an extensive 
study has been carried out also for several other decays [p0| . This can also 
be explained by a similar argument for the nuclei with Z or N = 40 — 50, in 
which the isolation of the Ogg/2 orbit plays an essential role. The 128 Te and 
130 Te cases [|19|, are a different matter. Another effect such as quadrupole 
deformation might be necessary to understand this case. 

The role of the shell structure has not been recognized sufficiently so 
far. Despite the arguments for the 48 Ca decay ||, there has been a lack of 
such discussions adapting to more general cases. In the 0vj3j3 decays, the 
difference of energy denominator among low- and high-lying intermediate 
states is much less important [0 than in 2vf3f3. However, the shell-structure 
mechanism remains and lower- lying GT strengths will still be significant, 
because of a certain similarity of the Qv(3f3 nuclear matrix-elements to the 
Ivfifi ones. 

The present calculation has been carried out with the algorithm developed 



in Ref.||10||. It has been confirmed extensively in this article how powerful 
this method is, by taking examples other than the 48 Ca decay as well as by 
a detailed discussion. There are a few other minor differences of the present 



study from that of Ref . |TIJ . While the GT operator proportional to at± has 
been used in Ref.|lOj, more realistic GT operators are adopted here. This 
does not influence the speed of convergence, as expected. Although only the 
computation starting from |GT~) is considered in Ref. JT0|, we have tested 



both the |GT+) and |GTj) cases. A slight faster convergence of the |GT^) 
case has become evident. It is naturally expected in the (3~ (3~ decays that 
the calculation starting from |GTj) converges somewhat faster than from 

|GTT>. 

The quasiparticle random-phase approximation (QRPA) has often been 
applied to calculating the f3j3 nuclear matrix-elements. It is known, however, 
that results of the 2vf3f3 calculation with QRPA are seriously sensitive to g pp 



around g pp = 1 [22) , where g pp denotes a factor introduced for the particle- 
particle interaction. This problem regarding g pp is crucial to predictability of 



the QRPA calculation. It has been shown in Ref.[P2]| that the low-lying GT 
spectra, particularly the GT + ones, are sensitive to g pp . Based on the above 
discussion of the shell-structure effect and the consequence of the present 
shell-model calculations, the g pp problem in the 2vf3(3 calculation seems to 
originate primarily in the sensitivity of the low-lying GT + strengths to g pp . 

It is emphasized that the low-lying GT strengths with E x < (5 — 10)MeV, 
rather than higher-lying component such as the GT~-resonance, are essential 
to the 2vf3f3 nuclear matrix-elements. There are several significant implica- 
tions. As is well-known, energies and wavefunctions of lower-lying states 
converge rapidly in the Lanczos diagonalization algorithm. Convergence of 
the present algorithm is accelerated further by the dominance of the low- 
lying GT strengths, because its essential part is the same as the Lanczos 
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method. The low-lying GT strengths will provide us with a crucial test on 
reliability of the calculations of the 2v(3(3 nuclear matrix-elements. It is true 
also for the QRPA calculations. Both of the parent and daughter nuclei of 
the (3(3 decays are stable against the single-/? decay. It is hopeful to measure 
by high-resolution charge-exchange experiments the low-lying GT strengths 
from the initial and final states. Theoretical approaches should be examined 
critically for the low-lying GT strengths by these measurements. As has been 
discussed in connection with Eq.(§), concise reproduction of individual GT 
strength is not necessarily required. It is requisite, however, to describe prop- 
erly the average character of the low-lying strengths in a sufficiently small 
energy range. This point is more important than to reproduce the total GT 
strengths. Moreover, if we had a reliable assessment on the signs of the GT 
matrix-elements, the 2v(3(3 matrix-element could be obtained experimentally 
by using the measured GT strengths. In the present calculations and those 
for the 48 Ca decay, energy regions adding positive and negative values to 
M w (2v(3(3) appear nearly separated. It may not be difficult to learn the signs 
from realistic calculations. 



4 Half-lives of 2v/3f3 decays 

The main purpose of this article is to clarify what is required to calculate 
the 2v(3(3 nuclear matrix-elements. As far as the calculated GT strengths 
from the initial and final states have not been examined sufficiently, we 
should not claim high precision for the 2v(3(3 matrix-elements calculated in 
the present study. We here show, however, the calculated half-lives, be- 
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cause there could be some interest in the prediction of the present realis- 
tic shell-model calculations. In order to evaluate half-lives of the 2vf3f3 de- 
cays, integration of the lepton degrees-of- freedom is carried out, by following 
Ref.g. As a result, half-lives of Tfy (EC/EC; 36 Ar -»■ 36 S) = 1.7 x 10 29 yr, 
Tff 2 (EC/EC; 54 Fe -> 54 Cr) = 1.5 x 10 27 yr are obtained from the present nu- 
clear matrix-elements. The half-lives of the 58 Ni decay is shown in Table [l], 
in comparison with the QRPA results^] and the data|23[]. The present cal- 
culation gives half-lives somewhat longer, but comparable to those by the 
QRPA calculation. Both calculations predict much longer lifetimes than the 
lower-limits of experiments. 



5 Summary 

It has been discussed what conditions are demanded to compute the 2v(3(3 
nuclear matrix-elements: We should prepare for a set of (approximate) inter- 
mediate states with (a) exhausting the relevant GT strength and (b) energy 
distribution having a sharp peak, over the energy eigenstates. An algo- 
rithm proposed in Ref.[[L(|] based on Whitehead's moment method seems 
quite promising from this viewpoint, particularly when it is associated with 
large-scale realistic shell-model wavef unctions. This algorithm has been in- 
spected, applied to the 36 Ar — > 36 S, 54 Fe — >■ 54 Cr and 58 Ni — > 58 Fe decays. 
Quite rapid convergence has been shown. It has also been pointed out that 
the low-lying part of the intermediate GT strength dominates the 2vj3j3 rates, 
owing to the shell structure of nuclei as well as to the energy denominator. 
However, it will not be always true that only the lowest intermediate 1 + state 
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is enough in evaluating the 2u(3(3 nuclear matrix-elements. It somewhat de- 
pends on the extent of the configuration mixing. Careful tests of the low-lying 
GT strengths, which will be doable with assistance of high-resolution charge- 
exchange experiments, are significant to predict the 2v(3(3 matrix-elements. 
Half-lives of those decays obtained within the present shell-model framework 
have been reported. 

The authors are grateful to Dr. S. Cohen for careful reading of the manuscript. 
The numerical calculation is financially supported by Research Center for 
Nuclear Physics, Osaka University (Project No.95-A-01). 
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Figure Captions 

Fig.|l]: Convergence of M w= q/ 2 (2z//3/3) for N. The upper sector shows the 
convergence starting from |GT^~) for the three decays. In the lower 
sector, the convergence starting from |GT+) and that from |GTj) are 
compared for the 36 Ar decay. 

Fig.|2j: GT + strength function for 54 Fe calculated by Whitehead's moment 
method: the strength functions obtained with N = 10 (upper), 30 
(middle) and 150 (lower) are compared. 

Fig.|3|: M w= q/ 2 (2^/3/3) with contribution of the intermediate states up to the 
excitation energy specified by the horizontal axis. 

Fig.^ GT + strength function of 36 Ar (upper), GT~ strength function of 36 S 
(lower) and M u)= q/ 2 (2uPP) up to E x for the 36 Ar decay (middle). 

Fig.||: Proton and neutron schematic configurations regarding the (3f3 decay. 
a,b): Ground-state configurations of the initial and final nuclei. 

c) Possible configurations generated by the GT + transition from the 
state of a). 

d) Possible configurations generated by the GT~ transition from the 
state of b). 



Tables 



Table 1: Half-lives (yr) of the 58 Ni — > 58 Fe decay. The values obtained in 
the present calculation are compared with those of the QRPA calculation of 
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